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ABSTRACT 



In efforts to demonstrate the linear Hubble law v = Hr from galaxy observations, the un- 
derlying simplicity is often obscured by complexities arising from magnitude-limited data. In 
' this paper we point out a simple but previously unremarked fact: that the shapes and orien- 

tations of structures in redshift space contain in themselves independent information about the 
cosmological redshift-distance relation. 

The orientations of voids in the CfA slice support the Hubble law, giving a redshift-distance 
power index p — 0.83±0.36 (void data from Slezak, de Lapparent, & Bijoui 1993) or p = 0.99±0.38 
(void data from Malik & Subramanian 1997). 

> 

Subject headings: large-scale structure of universe 

o . 

1. Introduction 

■ Hubble's (1929) observation that redshift increases linearly with distance for nearby galaxies has now 

been known for almost seventy years, and it is likely that its validity is not doubted by many. Yet, any attempt 



to demonstrate this simple law from galaxy observations soon descends into the complexities of magnitudc- 
limited observations, the broad galaxy luminosity function, Malmquist bias, and, at larger redshifts, K- 
corrections, evolution, etc. This has permitted an often-ignored but persistent challenge to a linear redshift 
law from those preferring a quadratic relation, with both challengers and supporters citing data ranging 
from relatively nearby, bright optical galaxies (Soneira 1979, Segal 1980) through the IRAS 1.2 Jy redshift 
catalog (Segal et al. 1993; Koranyi & Strauss 1997). 



One complication in understanding the classical Hubble diagram has been inhomogeneity, appearing 
as a dependence of density on radius or redshift or in clustering and peculiar velocities. With the advent 
of more galaxy redshift catalogs covering larger amounts of the sky, we obtain ever clearer pictures of the 
universe. The CfA "slice" (de Lapparent, Geller, & Huchra 1986), in particular, first revealed dramatic 
structures, large voids separated by well-defined walls, extending to a significant fraction of the survey 
volume, perhaps even calling into question whether a survey to this depth represents a statistically fair 
sample of the universe. In other ways, however, this inhomogeneity itself can be useful. In this paper we use 
the shapes and orientations of structures to obtain information about the background cosmology in which 
they are embedded. In Section 2 below we discuss how shapes and orientations of structures in redshift space 
depend on the redshift-distance relation, and in Section 3 we apply these considerations to the CfA slice (de 
Lapparent, Geller, & Huchra 1986) and the Las Campanas redshift survey (Shectman et al 1996, LCRS). 
Section 4 contains a final discussion. 
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2. Structure, Homogeneity, and the Redshift-Distance Relation 

The usual Hubble law cz — Hr relating cosmological redshift z and distance r is the linear case of a 
power-law relation characterized by an index p, 

z = Ar p . (1) 

(As in the usual relation derived from cosmological Robertson- Walker solutions to the Einstein equations, 
this behavior can be the small-r limit of a possibly more complex relation.) Since flux decreases with distance 
as r~ 2 , the classical Hubble diagram uses recession velocities and apparent magnitudes to infer p from the 
magnitude-redshift relation, 

5 z 

m = M + - log— . (2) 

V z 

However, this simple relation is difficult to observe in practice because of the broad distribution of galaxy 
luminosities. 

We present an alternative method making use of inhomogeneities to determine p from the shapes and 
orientations of structures in redshift space. Shapes and orientations in redshift space have been considered 
previously as indications of cosmological effects that depend on the density parameter f2 , the deceleration 
parameter qo, or the cosmological constant A, using as tracers the distribution of galaxies in clusters (Alcock 
& Paczynski 1979), the angle with the line of sight of the line joining neighboring pairs of quasars (Phillipps 
1994), or distortions of voids (Ryden 1995). Because these effects are weak at small redshift, with typical 
behaviors of the form 1 + kz for some k, implications for cosmological parameters have been inconclusive. 
However, the effect of a redshift index p ^ 1 appears at zero order. 

For example, suppose that we can identify features with a characteristic angular extent AO and redshift 
extent Az. From the redshift-magnitude relation of equation ([!]), A8 is related to the transverse physical 
diameter of the structure r± by AO = Ar±/r; while Az is related to the line-of-sight diameter Arn by 
Az/ z = pArn/r (Ar <C r). Thus, without ever needing to know the distance r, we have 

1 Az Ar,, 

~zA9 =P A^> (3) 

the aspect ratio in redshift space is related to the aspect ratio of the physical object by exactly the redshift 
index p. Ignoring peculiar velocities, a (distant) sphere is mapped in redshift space into an ellipse stretched 
by a factor p along the line of sight. 

With irregularly shaped objects, sparsely sampled, this aspect ratio may be difficult to measure accu- 
rately. Another quantitative measure of p can be obtained from considering orientations. The line joining 
the center of a sphere to a point on the surface makes an angle with the line of sight (this is always an 
angle between and 90°). Consider the average of the cosine of this angle in redshift space, = (cos 6*). In 
three dimensions this is simply 

M3d - -4t- (4) 
p+1 

Inverting, the redshift index p is found from /i as p = /i/(l — /x), with error Ap = 0^/(1 — /i) 2 . For a 
two-dimensional sample such as the CfA slice, the expression is a bit more complicated, 

2 / ^2 

fJ>2d=-\ -n -tan -1 vV - 1 (5) 

7T V P — 1 
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for p > 1, analytically continued for p < 1, with /^2<2 = 2/tt and dfi/dp = 2/3tt at p = 1. For all p we have 
M2d > M3d, but the basic behavior remains the same, /i is a smooth, monotonic function of p. Figure 1 shows 
\i as a function of p for both two and three dimensions. These expressions hold not only for physical objects 
that are inherently round, but also for ellipses of any eccentricity, provided only that orientations to the line 
of sight are random. Equation rtq) will allow us to determine p from data. 



3. Application to Observations 

We first present the effect visually using the Las Campanas Redshift Survey (Shectman et al. 1996), 
comprising in total 26,418 redshifts in six 1?5 x 80° strips in selected north and south pole areas with 
median redshift of approximately 30, 000 kms -1 . Figure 2a shows the LCRS —12° declination slice in angle- 
redshift space, which is also physical 6-r space if distance and redshift are related by the linear Hubble law. 
Represented in the figure are 4071 LCRS galaxies with cz less than 45, 000 km s -1 . Figure 2b shows how the 
slice would look in 6-r if r ~ z 2 (p = A in eq. B; this is also how the slice would look in redshift space if 
true distances are as in Figure 2a and p = 2). Figure 2c shows physical distances if r ~ z 1 / 2 (p — 2; this is 
how the slice would look in redshift space if distances are as in Figure 2a and p = Aj. In Figure 2a we see 
some structures that are elongated radially, and some that are elongated in the transverse direction, as we 
expect if the universe is homogeneous and we do not live in a privileged location. In Figure 2b, however, 
essentially all elongations are radial, and in Figure 2c, essentially all are transverse. Further, only for p = 1 
is the density not steeply varying with r. 

We can attempt to put this on a quantitative basis by examining nearest neighbor orientations in the CfA 
"slice of the universe" (de Lapparent, Geller, & Huchra 1986), 1093 galaxies brighter than 15.5 in magnitude 
in the 6° declination band 25? 5 < S < 32? 5 covering 135° in right ascension from 8 h to 17 h . For the 1057 
galaxies in the inner 80% of the CfA slice, the average of the orientation of the line joining nearest neighbors 
gives fi = 0.644 ± 0.010. From eq. (||), this corresponds to p = 1.037 ± 0.050. This appears to support the 
Hubble law, but in fact is only a reflection of the metric we have used to determine "nearest neighbor" in 
redshift space, As 2 = Az 2 + z 2 A6> 2 . If we first square all the redshifts and repeat the analysis we obtain 
fi = 0.670 ± 0.010; if we take the square root of all redshifts we obtain fi = 0.634 ± 0.010. The explanation is 
simple: the radial rescaling changes which points are nearest neighbors in redshift space. This may at first 
be surprising, but seems obvious after considering the following. Take a uniform distribution of points in the 
x-y plane and then stretch it by an arbitrary amount along the x-axis. The resulting distribution was, and 
thus still is, independent of y; and it was, and so still is, independent of x; and although the overall density 
may be changed, the distribution must again be uniform in the portion of the plane it occupies. Thus we 
conclude that nearest neighbors do not provide the quantitative measure of p we seek. 

We need to use objects whose identity is not changed by remapping the radial coordinate, and an obvious 
candidate is voids. The identification of a void relies only on distinguishing locally underdense regions and not 
on an absolutely calibrated sensitivity. Isolated voids become innately round; unlike collapsing objects, which 
typically collapse along one axis first, non-spherical voids become more spherical as they evolve (Bcrtschinger 
1985; Ryden 1994). Even if a typical bubble is not perfectly round in shape, from homogeneity we expect 
the distribution of shapes to be independent of z, and from isotropy we expect orientations to be random. 
We use the orientations of voids identified in two wavelet analyses of the CfA slice to determine the redshift- 
distance index p. Slezak, de Lapparent, & Bijoui (1993) present a list of 15 voids, identified using an 
isotropic template with orientation defined by fitting to an ellipse after their wavelet analysis detects a 
significantly underdense region. Malik & Subramanian (1997) list a total of 29 voids, from which we remove 
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three apparent duplications and one that is perfectly round, so that no orientation is determined. Their 
anisotropic wavelet templates includes orientation as a parameter, but their angular resolution is limited to 
multiples of 30°. Results are given in Tables 1 and 2. The best determinations are those averaged over all 
15 Slezak voids, /i = 0.597 ± 0.0941 corresponding to p — 0.831 ± 0.357; and the la significance level from 
Malik & Subramanian, which gives [i = 0.635 ± 0.0807, or p = 0.992 ± 0.377. Values for various subsets are 
given in Table 1 and Table 2. If we look for the best integer index, these results prefer p = 1 to p = 2. 

We might expect peculiar velocities to distort the apparent shapes of expanding voids. Voids generally 
expand more rapidly than the Hubble flow, but the expansion velocity of the shell of material around a hole 
typically is only a small fraction of the Hubble velocity across the hole (Peebles 1982; Ryden 1994). For 
isolated, spherical "top-hat" voids we can compute the expansion velocity, typically expressed as a fraction F 
of the Hubble flow. For ft = 1, this factor is F = ^ if the void is uncompensated and F = ^ if compensated 
(Bcrtschinger 1985), which becomes F = ifi ' 6 for f2 <; 1 (Regos & Geller 1991); however, the flow at a 
void boundary generally depends critically on its surroundings. Measured differences between redshift and 
independent Tully-Fisher distance determinations for galaxies in the largest bubble of the CfA slice confirm 
that peculiar velocities are small, consistent with zero and with a Ter upper bound less than 5% of the 
diameter (Bothun et al. 1992; this determination of course assumes p = 1). It is instructive to compare 
with the results for overdense regions, also provided by Malik & Subramanian (1997). As a whole, overdense 
detections are smaller and more anisotropic than are the underdense regions, and for the higher significance 
level detections, they are preferentially oriented along the line of sight, with /i = 0.827 ±0.066 for significance 
> lOcr. This is as we expect; collapsed objects are smaller and have larger internal velocity dispersion. 



4. Discussion 

We have shown in this paper how the apparent shapes and orientations of objects in redshift space 
can be used to determine the redshift-distance relation. When redshift is not linear in distance, objects 
such as voids that are intrinsically round in space when viewed in redshift space are stretched along the 
line of sight; and an initially isotropic distribution of orientations becomes distorted in the radial direction, 
with measurable effect, as in equation (|^). Application to voids in the CfA slice gives redshift power index 
p = 0.83 ± 0.36 from Slezak et al. (1993) and p = 0.99 ± 0.38 from Malik & Subramanian (1997), both a 
modest preference for p = 1 over p = 2. To obtain these results in redshift space with a true p = 2, voids in 
the CfA slice in space would have to be flattened and preferentially aligned transverse to the line of sight. 
The main complication to the simple interpretation is likely to come from peculiar velocities; an expansion 
velocity will introduce an apparent line-of-sight elongation. However, this effect is expected to be small: 
in the most extreme case, of an uncompensated, completely empty, isolated void in an f^o = 1 universe, 
expansion would make p — 1 appear to be p = |. In any case, unless voids are contracting, a condition 
difficult to make physical sense of, peculiar velocities will only increase, not decrease, the apparent value of 
p, and will not distort a quadratic redshift law to appear to prefer p = 1 . 

We do not expect void orientation statistics to replace classical methods for demonstrating the Hubble 
law. The most precise confirmation of the linearity of the Hubble law, p/5 — 0.2010 ± 0.0035, has been 
obtained recently from the Hubble diagram of type la supernovae (Riess, Press, & Kirshner 1996). The 
luminosity function as a function of redshift and the radial dependence of galaxy density in the IRAS 1.2- 
Jy catalog (Koranyi & Strauss 1997) also support the linear Hubble law. However, our method offers an 
independent verification of the Hubble law in which magnitude-limited data, Malmquist bias, and evolution 
do not present problems. Further examination of the effect in real data, especially including in detail the 
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peculiar velocities of the void walls, will undoubtedly introduce complications to the simple relations in 
equations (Q) and (||). Whatever these complications, they will be different from those that enter arguments 
over the redshift-magnitude relation. 

We thank MJG for providing the CfA slice data in machine readable form. This research was supported 
in part by NASA grant NAG5-2835. 
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Tabic 1. Rcdshift Index p from Slezak et al. (1993) Voids 



Void size (kms 1 ) 


N 




P 


350-700 


7 


0.594 ±0.130 


0.821 ±0.488 


1000-2000 


8 


0.598 ±0.143 


0.836 ±0.547 


350-2000 


15 


0.597 ±0.094 


0.831 ±0.357 



Table 2. Rcdshift Index p from Malik & Subramanian (1997) Voids 



Significance 


N 




P 


> 6cr 


22 


0.685 ±0.070 


1.265 ±0.477 


> 7(7 


18 


0.635 ±0.081 


0.992 ±0.377 


> 8cr 


8 


0.654 ±0.122 


1.086 ±0.634 
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Fig. 1. — The average /i = {cos 9) over orientation relative to line of sight in redshift space for an isotropic 
distribution of orientations in space, as a function of the redshift index p. The dashed line shows the result 
in three dimensions, the solid line the result for a two-dimensional slice. 




Fig. 2a. — The LCRS galaxy distribution in z and 9. 




Fig. 2b. — The LCRS galaxy distribution in r and 9 for z ~ r 1 / 2 , or r ~ z 2 . 




Fig. 2c— 



The LCRS galaxy distribution in r and for z ~ r 2 , or r ~ z 1 / 2 . 



